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A hexagonal warping term has been proposed recently to explain the experimentally observed 2D equal 
energy contours of the surface states of the topological insulator Bi2Te3. Differing from the Dirac fermion 
Hamiltonian, the hexagonal warping term leads to the opening up of a band gap by an in-plane 
magnetization. We study the transmission between two Bi2Te3 segments subjected to different in-plane 
magnetizations and potentials. The opening up of a bandgap, and the accompanying displacement and 
distortion of the constant energy surfaces from their usual circular shapes by the in-plane magnetizations, 
modify the transverse momentum overlap between the two Bi2Te3 segments, and strongly modulate the 
transmission profile. The strong dependence of the TI surface state transport of Bi2Te3 on the magnetization 
orientation of an adjacent ferromagnetic layer may potentially be utilized in, e.g., a memory readout 
application. 



Topological insulators (TIs)^ are a relatively new class of materials which have attracted much attention due to 
their rich physics. For example, Majorana fermions have been shown to exist in junctions between super- 
conductors and ferromagnetic insulators deposited on top of TIs^"^. Furthermore, strong spin orbit coupling 
leads to interesting effects when a magnetization is applied to a TP"^. In particular, Yokoyama and colleagues 
studied the magnetoresistance in a two dimensional junction between two ferromagnets magnetized in different 
directions deposited on top of a topological insulator^ They found that the transmission at a given energy is 
strongly influenced by the relative /c- space displacements of the Fermi surfaces resulting from the differing 
magnetizations. The influence of the Fermi surfaces on the magnetoresistance motivates the study of the topo- 
logical insulator Bi2Te3 in this present work. 

The experimentally observed^'^° 2D equal energy contours (EECs) of the surface states in Bi2Te3 differ at high 
Fermi energies from the circular Dirac cone predicted by the simple Dirac fermion Hamiltonian H = vkx a. At 
low Fermi energies, the EEC takes the form of a circle (Fig. 1(a)). As the Fermi energy increases, the contour 
evolves from a circle to a hexagon and then to a snowflake with sharp tips along the six TM directions. Based on 
the underlying 3-fold rotational and two-fold mirror symmetry of the [111] surface of the underlying rhom- 
bohedral Bi2Te3 crystal structure, Fu suggested the addition of a hexagonal warping term to the Hamiltonian^ \ 
The Hamiltonian then takes the form of 



H = v(kx<Jy 



(1) 



where x is in the YK direction. The Hamiltonian reproduces the experimentally measured EECs. 

Unlike the Dirac cone which is rotationally symmetric about the k^, ky plane, the dispersion relation arising 
from Eq. 1 has only 6 fold rotational symmetry. The anisotropy of the Fermi surface in the k^, ky plane can be 
expected to give rise to interesting directional dependent effects. For example, it has been suggested that a Bi2Te3 
segment with a potential step may be exploited as a flat lens for focusing electron beams^^. 

Another peculiarity of the Hamiltonian is that an in-plane magnetic field modeled by a Zeeman term can open 
up a band gap^\ The Zeeman term mathematically resembles that due to magnetic coupling to a proximate 
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Figure 1 | The EECs and eigenstate spin orientations for E = (a) 0.1 eV, (b) 0.3 eV and (c) 0.5 eV. The green arrows on the EECs indicate the 
orientation of the (a^ and (cr^) components at each point on the contour with larger arrowheads indicating higher spin polarization. The red and blue 
dots indicate the (cr^) polarization with larger dots indicating larger magnitudes of spin polarization. Blue (red) dots indicate positive (negative) spin z 
polarization. The dotted line in panel (c) indicates one value of ky where there are six real values of 



ferromagnetic film or the effects of magnetic doping^^. For ease of 
exposition we say that the Zeeman term arises from an external 
magnetization. We write the Zeeman term as M-a and collect the 
appropriate coupling constants into M itself. The introduction of a 
bandgap due to an in-plane magnetization differs from the simple 
Dirac fermion Hamiltonian H = vkxd where only an out-of-plane 
magnetization can open up a band gap. We find that besides opening 
up a band gap and displacing the position of the Dirac points in 
reciprocal space, the in-plane magnetic field also distorts the shape 
of the constant energy surfaces. 

Despite the interesting bandstructure properties of Bi2Te3, trans- 
port across Bi2Te3 segments has not yet been studied extensively. 
Ref. 14 derived the transmission between one TI segment and 
another with an externally applied step -edge potential V(x) = 
— Vo6(—x) where 0(x) is the Heaviside step function. In the deriva- 
tion it was implicitly assumed that there is a propagating state in the 
transmitted segment which satisfies energy and transverse 
momentum conservation. However, there exists more than one 
propagating state with real, positive for some values of ky at high 
\E — V\ (see Fig. 1(c)). It is not evident which of these multiple values 
of /c^c in the transmitted region should be chosen. Moreover, the 
derivation matched only the wavefunctions, and not the probability 
flux, across the interface. Ref. 16 considered the transmission of a 
wave incident from x<0 across a Dirac delta potential located at x = 
0. The transmitted wavefunction at x > 0 contains left propagating 
components back into the incident region within the sinik^x) and 
cos(kxX) terms in Eq. (4) of the paper. The presence of backwards 
propagating wave components in the transmitted region is counter- 
intuitive. 

For a given value of transverse momentum ky and Fermi energy E, 
the Hamiltonian Eq. 1 admits six eigenstates with real or imaginary 
kx which may be real or complex, and which come in three Kramers 
time-reversal pairs. We consider the transmission of electrons travel- 
ing in the +x direction from one Bi2Te3 segment on the left to 
another on its right. Three of the eigenstates propagate or decay in 
the +x direction, and are physically appropriate states in the segment 
which a wave traveling in the +x direction is transmitted into. The 
other three eigenstates propagate or decay in the —x direction, and 
are physically appropriate as the reflected components of a wave- 
function in the incident segment. Both Refs. 14 and 16 assume that an 
incident energy eigenstate with wavevector k^ be reflected at the 
interface into an eigenstate with the wavevector —k^. This imposes 



an artificial constraint that the incident wave is not reflected back 
into the other 2 eigenstates which also propagate or decay in the 
correct direction. 

One difficulty in considering all eigenstates in studying transmis- 
sion across an interface between 2 Bi2Te3 TI segments is the issue of 
boundary conditions. In each TI segment there are 6 eigenstates for a 
given value of E and transverse momentum ky. In solving the trans- 
mission problem across an interface from TI region I to region II, the 
weightages of six of the twelve eigenstates involved are fixed by the 
requirements that evanescent states decay away from the interface, 
and that the state incident on the interface from region I is given. This 
leaves the weightages of the remaining 6 states to be solved for. The 
usual practice of matching the wavefunctions and first derivatives for 
the two spinor components across an interface between the two TI 
segments gives only 4 equations. Another boundary condition 
besides wavefunction and first derivative continuity is required to 
solve for the 6 unknowns uniquely. We note that more recent 
works^^'^^ did consider all six eigenstates explicitly. However, Ref. 
15 only considered transmission across a delta potential barrier, 
whereas it is not clear from the paper what set of boundary condi- 
tions Ref. 17 used. 

To the best of our knowledge, transport across Bi2Te3 segments 
magnetized in different directions has not yet been studied. In this 
work, we show that probability flux conservation across the interface 
requires that the matching of the second derivatives of the wavefunc- 
tions across the interface. We then calculate the transmission 
between two Bi2Te3 segments subjected to different external poten- 
tials and magnetic fields. We find that the distortion and displace- 
ment of the constant energy surfaces, and the opening up of a band 
gap lead to a rich transmission profile as the magnetization directions 
and energy are varied. 



Methods 

Boundary conditions. We consider a generic Hamiltonian of the form 

= - iad^ - bdl + icdl + d. 



(2) 



The a, b, c and d above may be numbers or Hermitian operators which commute with 
px, for example some combinations of ? or py and are independent ofp^. For this 
dH 

Hamiltonian, we have Vx = = a + 2bpx + 3cp^. Explicitly expanding out dtp gives, 
after some manipulations of the x derivatives and integration by parts. 
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Figure 2 | (a) The band gap in eVas a function of the magnetic field strength Mand angle 0m- (b) The thick blue lines indicate the trajectory of 
the Dirac points in k space as 0m is varied at fixed values of Mfrom 0.05 eV to 0.20 eV. The black dotted lines show the constant energy surfaces at 
E= 0.1 eVand£= 0.2 eV. 



dtP = i[{HiJyyily-r{HiJy)) 

= -d^ {ij/*ail/ + 2Im(0* bd^il/) - 2Re(0*ca^0) + {d^il/*)c{d^il/)) . 



(3) 



Comparing with the continuity equation dtp + d^j = 0, we recognize the terms within 
the large brackets in the second line as the probability flux. If c = 0 (i.e. the 
Hamiltonian is only quadratic inpx), the flux is the real part of \J/*v\J/ so that for a 
Hamiltonian quadratic inpx, matching \J/*v\J/ across an interface - which matches 
both the real and imaginary parts - automatically matches the flux across both sides of 
an interface. 

However, for ct^O, the real part of the c term in if/* vij/ is 3Re (xj/* cdlij/) which is not, 
in general, equal to the c term in the flux, 2Re(i/^*c5^i//) — dxijy*cdxijy. A simple 
example suffices to illustrate this. Putting c = 1 for simplicity, and writing \J/ = 
expiikx) where k = {kr + iki) is in general complex, we have at % = 0, we have 
3Re{iJ/*cdliJ/) = -?>(k]-k]) which is not equal to 2Re(i/^*ca^i/^) -dxiJ/*cdxiJ/ = 
— 3k^ + 2k^ if the wavevector has an imaginary part. 

In the system we will be considering in this paper, we identify the operators 
a = VGy — ?)k^ylGzy b = 0, and c = Xo^ across both sides of the interface. (They differ in 
the d term of Eq. 2.) In this instance, flux conservation can be achieved by matching i) 
x\i ii) dxii and iii) d\x\i across both sides of the interface. 

Results 

Equal energy contours under an in-plane magnetization. The 

eigenenergy E for a given ^ = A;(cos((/>),sin((/>),0) and magnetiza- 
tion M=M(cos((/>M)?sin((/>^),0) reads 

£2 _ i (A;^^^ 1 + cos (6(/> ) ) + 2/:^ + IM^ - 4M kv sin((/> -^J)) .(4) 

The high powers of k and presence of trigonometric terms involving 
^ and present difficulties for the derivation of compact analytic 
expressions for the band gap and position of the Dirac point in k- 
space at which the lowest energy propagating particle state with real k 
exist. We calculate these quantities numerically. We use the 
numerical values of v = 25.5 eVnm and X = 0.25 eVnm"^ ^\ 

The in-plane magnetization leads to the opening up of a band gap 
and a displacement of the Dirac point, as shown in panels (a) and (b) 
of Fig. 2 respectively. The band gap reaches its maximum value of M 
at magnetization angles (/>m = ±{2n — l)7i/6, n= 1, 2, 3, and varies in 



a roughly sinusoidal manner with the magnetization angle. The band 
gap shifts the minimum energy at which propagating particle states 
exist upwards. This leads to the /c-space area covered by the EEC at a 
given energy being generally larger for a magnetization angle with a 
smaller band gap compared to that of a magnetization angle with a 
larger band gap. This fact will have consequences for charge trans- 
port between Bi2Te3 segments magnetized in different directions. 

The trajectory of the Dirac points for a given magnetization mag- 
nitude as the magnetization angle is varied follows a more compli- 
cated path than the circular trajectories of the simple Dirac fermion 
Hamiltonian. For the latter, the trajectories form circles in A;-space. In 
Bi2Te3, the trajectories of the Dirac points resemble the shapes of the 
EECs. The trajectories of the Dirac points exhibit greater distortion 
from the circular shape at a given value of M than the EECs do for 
E = M. 

The positions of the Dirac points k* differ from the k* = zxM/v 
given in Ref. 11. The latter holds only in the low M limit when the 
trajectories of the Dirac points are roughly circular. The /c-space 
angle (/> = eiYctan{ky/kx) of the Dirac point is related to the magnet- 
ization angle by (/>m — ^ ~ 7^/2. This can be seen from the last term of 
Eq. 4 where the energy achieves its minimum value for fixed values of 
k and M when the sin((/> — (/>m) term is 1. At (/> = (/>m + 7r/2 and the 
large k limit, the k^ term of Eq. 4 dominates so that E ~ k^X 
|sin(3(/>M)|- This explains the peaking of sinusoidal variation and 
peaking of the band gap at (^m — — (2w — l)n/6, n = l,2,3. 

Transport between two Bi2Te3 segments. To study the effects of the 
in-plane magnetization on the transport properties of Bi2Te3, we 
consider the system outlined schematically in Fig. 3. The system 
consists of two Bi2Te3 segments of semi- infinite length along the x 
direction and infinite width along the y direction. Charge flows from 
the left source segment to the right drain segment in which a 

magnetization M2 =M2(cos((^^2)'Sin((^^2)) ^ potential U is 
applied. The potential shifts the dispersion relation of the drain 
segment on the energy axis with respect to that of the source 
segment. For a given source Fermi energy E we define E2 = E — U. 
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Figure 3 | A schematic of the system studied. The system consists of two semi-infinite long Bi2Te3 segments of infinite width. Charge flows from the left 
source segment into the right drain segment in which a magnetization M2 and external potential 17 is applied. 



Fig. 4 shows the transmission plotted as a function of magnitude 
M2 and direction (^^2 of the drain magnetization for the transmission 
from an unmagnetized source segment. The electron energy is cho- 
sen to be the intermediate range, i.e., E = 0.3 eV, at which the EEC 
takes the form of a hexagon, while the drain segment is biased at E2 = 
E — U2 = 0.1 eV, i.e., at the low energy regime where the EEC 
assumes a circular shape. 

The dependence of the transmission on M2 can be divided into 
three regimes. In the low M2 regime (0 eV < M2 < 0.07 eV) the 
transmission peaks weakly around (/>m2 = —nil, while at intermedi- 
ate values of M2 the transmission exhibits sharper peaks near (/>m2 — 
0, ± 71. At large values of M2, the trend is reversed and the transmis- 
sion drops to almost zero near (/>m2 ^ ±n7i/6, n = (1, 3, 5). 

Discussion 

These observations may be explained by examining the evolution of 
the EECs as M2 and (/>m2 are varied (as shown in Fig. 5). 

For low values of M2, the EECs at the drain lie completely within 
the /c-space region spanned by the EEC at the source. Fig. 5(a) shows 
that the drain EECs adopt an almost circular profile and the ky range 
spanned is only weakly dependent on (/>m2- The spin polarizations 
(not shown) of the right-propagating eigenstates also do not exhibit 
much angular variation. Thus, there is little modulation of the trans- 
mission as the orientation of M2 is varied. In the intermediate range 



of M2 [shown in Fig. 5(b)], the transmission is dominated by the ky 
range spanned by the drain EEC. Since the transverse momentum is 
conserved across the source-drain interface, a larger range of ky 
subtended by the drain EEC translates to a greater number of 
propagating states which can transmit electrons over a wider range 
of incidence angle. The transmission peaks occur near (^^2 — 0, n, 
which correspond to the two magnetization angles where the drain 
EECs cover the largest ky range. In the large M2 regime shown in 
Fig. 5(c), the energy gap induced by the in-plane magnetization 
eliminates the presence of propagating states at (/>m2 — — nn/6, n = 
(1, 3, 5), and reduces the ky range covered by the drain EECs at other 
magnetization angles. The absence of propagating states for certain 
range of (/>m2 leads to a sharp modulation of the transmission about 
the critical angles of (^m2 = ±n7i/6, n = (1, 3, 5), a feature which 
would be conducive for applications to sense the orientation of M2, 
such as memory read-out. The displacement of the drain EECs from 
the /c-space origin also leads to some portions of the drain EECs lying 
outside the source EEC. This leads to the generally lower transmis- 
sion seen in Fig. 5 for large values of M2. 

The effects of varying the potential in the drain segment for a fixed 
magnitude of magnetization provides an alternative perspective. 
Fig. 6 shows the transmission plotted against the drain magnetization 
direction and drain potential. The trends observed can also be 
explained by examining the EEC profiles in the drain region, shown 
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Figure 4 | The transmission plotted against the drain magnetization direction ^mi and magnetization magnitude for E = 0.3 eV and £2 = E— 1/2 = 
0.1 eV. 
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Figure 5 | The EECs in the source region (thick black hexagon) at £ = 0.3 eV, the trajectories of the Dirac points traced out by by the variation of 
^M2 (thick green lines) and the EECs in the drain region (thin blue lines) for E2 = 0.1 eV and M2=(a)0.05 eV, (b) 0.15 eV and (c) 0.28 eV. The 

positions of the Dirac points and the corresponding values of for each of the drain EECs are marked out in (b) in the same color. The ky range spanned 
by the = 0 drain EEC is indicated in panel (b). 



for three representative values of E2 in Fig. 7. The low transmission 
centered around (/>m2 = mil6, n = 1, 3, 5 at low values of £2 corre- 
sponds to the absence of propagating states due to the large band 
gaps induced at these values of (/>m (Fig- 7(a)). The transmission peaks 
at = 0, 71 at intermediate values of where the drain EECs still lie 
largely within the ky region spanned by the source EEC, correspond 
to the largest range of ky spanned by the drain EECs as 0M2 IS 
varied (Fig. 7(b)). As £2 is increased further, an increasing proportion 
of the /c-space area spanned by the EECs starts to fall outside the ky 
region spanned by the source EEC. This leads to the transmission 
peaks occurring close to (^m — ±7i/2 where the overlap in the ky 
region spanned by both the source and drain EECs is the greatest 
(Fig. 7(c)). 

We have thus far concentrated on the ky overlap between the 
source and drain EECs as the primary explanation of the observed 



transmission profile. However, the k^ and spin mismatch in the 
source and drain regions does affect the transmission as well. In order 
to investigate the effects of these two factors, we now let the range of 
E2 in Fig. 6 extend to negative values as shown in Fig. 8. Negative 
values of £2 correspond to transmission from source particle states to 
drain hole states, which have the opposite spin configuration. 

The general features of the transmission profile for particle-to- 
particle transmission (positive £2) are largely similar, but not ident- 
ical, to those for particle-to-hole transmission. The reason for this Ues 
in the fact that for a given magnetization M and |£| the EECs for the 
particle and hole states with energy ± |£| have identical shapes. The 
right propagating states for positive and negative values of £, how- 
ever, lie in different segments of the EECs as shown in panel (b) of 
Fig. 8. The ky range spanned by the right propagating states for 
negative and positive values of £ are roughly the same. This accounts 
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Figure 6 | The transmission plotted against the drain magnetization direction direction ^m2 and drain energy £2 for £ = 0.3 eV and M2 = 0.2 eV. 
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Figure 7 | The EECs in the source region (thick black hexagon) at £ = 0.3 eV, the trajectories of the Dirac points traced out by the Dirac points 
(thick green lines) by the variation of (green dots) and the EECs (thin blue lines) in the drain region for M2 = 0.1 eV, and £2=(a)0.05 eV, 
(b) 0.12 eV and (c) 0.25 eV. 
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Figure 8 | (a) The transmission profile for the same set of parameters as in Fig. 6 with the exception that £2 now extends to negative values, (b) The 
EECs for an unmagnetized segment at £ = 0.1 and no magnetization (black), and M2 = 0.2 eV, — 0-47i and 0m2 — ~0.47i respectively as 
labeled. The arrows indicate the directions of the in-plane spin polarization on the EECs of the right propagating hole (green) and particle (red) states in 
the drain segment. The lengths of the arrowheads are indicative of the magnitude of the in-plane spin polarization. 
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for the similarities in the particle-to-particle and particle-to-hole 
transmission profiles in terms of the rough positions of the trans- 
mission peaks and troughs. However, the difference in the kx values 
and spin orientations of the right propagating states of electrons and 
holes leads to the asymmetry in the particle-to-particle and particle- 
to-hole transmission profiles. This asymmetry is accentuated by the 
distortion of the constant energy surfaces from the circular shape by 
the hexagonal warping term. In addition, there is also an asymmetry 
in the transmission profile about the (^m2 — 0 axis. This may be 
explained by considering the EECs for ±(/>mj for some given (/>m- 
As shown in Fig. 8(b), the EECs for ±(^m are reflections of each other 
about the = 0 axis. However, the portions of the EECs which 
correspond to the right propagating states lie on different regions 
of the corresponding EECs for ±(^m- This leads to, for instance, the 
right propagating states at a given value of ky having different spin 
polarizations for opposite signs of (/>mj and hence an asymmetrical 
transmission profile about (^m2 = 0- Thus, in summary, two factors 
play a role in determining the dependence of the transport on the 
magnetization direction and magnitude - the extent of the overlap 
across ky of the the source and drain EECs, and the degree of mis- 
match between the spin configurations of the propagating states in 
the two electrodes. 

We note in passing that proposals for the formation of the 
Majorana fermion state in TI systems involve junctions between 
ferromagnetic insulators and superconductors deposited on top of 
TIs. The anisotropy in the shape of the EECs induced by the in-plane 
magnetization and asymmetry between the particle and hole states 
just discussed may lead to the presence of Majorana fermion states in 
Bi2Te3 with unusual properties compared to Majorana states result- 
ing from TIs with the more commonly studied circular EECs. The 
anisotropy may also be indicative of anomalous effects which emerge 
as the dimensions of the TI segments are shrunk from being semi- 
infmite (in the x and z directions) to finite ones^^. A detailed discus- 
sion of these issues are, however, beyond the scope of this paper. 

Conclusion 

In this work we showed the appearance of a band gap, and the 
accompanying distortion and displacement of the constant energy 
curves in k space in Bi2Te3 with the application of an in-plane mag- 
netic field. We then studied the transmission to a drain Bi2Te3 seg- 
ment with an external magnetization and potential applied from a 
source Bi2Te3 segment without the magnetization and potential. The 
band gap, and distortion and displacement of the constant energy 
curves by the in-plane magnetization affects the overlap of the trans- 
verse momentum ky range spanned by the EECs in the source and 
drain regions. The sharp modulation of the transmission with the 
magnetization orientation for the high M2 regime and low E2 regime, 
as depicted in Figs. 4 and 6, respectively, suggests a possible applica- 
tion of Bi2Te3 in the read-out of magnetic memory. The two states of 
a memory bit can be represented by the magnetization direction of an 
adjoining ferromagnetic layer (which is coupled to the topological 
surface states of Bi2Te3) being in the (/>m ^ 0 and (/>m = 7i/2 directions. 
In the former state there is no band gap and the transmission is finite 
while in the latter the transmission is 0 due to the large band gap. A 
gate voltage can be applied on the drain segment to tune the value of 
E2 so as to achieve an optimum balance between the differences in 
transmission between the two memory states, and the robustness of 



the system to deviation of the magnetization from the reference (/>m ' 
0 and — 7i/2 directions. 
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